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SUMMARY 

The recently derived particular integrals of the Prandtl-Busemann 
iteration equations mak e possible the extension of the familiar source- 
sink concept to the solution of the higher-order iteration equations for 
the subsonic potential, flow over thin sharp-nose symmetric two- 
dimensional profiles. An explicit expression is derived for the second- 
order velocity potential and velocity components' and a method for 
obtaining the higher-order terms is indicated. The velocity at the sur- 
face of the Kaplan bump is evaluated to illustrate the method. 


INTRODUCTION 


Iteration methods have been used extensively during recent years 
for the calculation of compressible flows. The Prandtl-Busemann or 
Ackeret iteration process is the one most often applied in aeronautics 
and is the one considered herein. In the past this iteration technique 
has been used primarily to calculate the flow past specific profiles, 
each profile presenting a distinct problem involving a great amount of 
labor in its solution. 

Van Dyke (reference 1) has given the particular integral for the 
second-order Prandtl-Busemann iteration equation and recently Kaplan- 
(reference 2) has obtained the particular integrals for both the second- 
order and third-order equations and has described a method for obtaining 
higher-order ones as well. With the use of these particular integrals, 
the higher -order solutions for the flow over sharp-nose symmetric bodies 
in two-dimensional flow are easily obtained by the method to be described 
herein, and the calculation of the surface velocity or pressure distri- 
bution for a specific profile requires only an integration. The procedure 
to be described involves considerably less labor than formerly required 
and also has the advantage of simplicity of treatment since the familiar 
concept of a source-sink distribution in the plane of symmetry is 
employed in the analysis. 



2 


MCA TN 2253 


A general expression for the second-order velocity potential for 
symmetric profiles is developed and, as an example, the flow over the 
Kaplan hump is treated in detail. The procedure may be extended to 
obtain the third- and higher-order perturbation velocities. 


ANALYSIS 


The Prandtl-Busemann iteration equations for obtaining approximate 
solutions to the exact nonlinear equation governing the two-dimensional 
flow of a compressible fluid are briefly developed as follows 
(reference 2): 

The differential equation for the velocity potential $ is 
( X " ^2)^ " 2 ^2 *xy + ^ - ^yy = 0 

where 

* 

x,y rectangular Cartesian coordinates in flow plane 

u, v velocity components along x- and y-axes, respectively 

a local speed of sound 

and the subscripts x and y denote partial differentiation with 
respect to these designated variables. With the introduction of a 
characteristic length c/2, where c. is the chord, and the undisturbed 
stream velocity U as the unit of velocity, the quantities x, y, u, 

v, and <I> are used hereinafter to denote, respectively, the nondimen- 

sional quantities 2x/c, 2y/c, u/u, v/u, and 2$>/Uc while a retains 

its original meaning. 

Let 

Y = tY x (x) + t^gU) + 0(t3) , (-1< x < l) 

and 

Y = 0 (x > lj x $ -1) 


where 0(t3) denotes terms of the order t^, define a slender sym- 
metric body of thickness coefficient t lying on the x-axis between -1 
and 1. 
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The appropriate boundary conditions for the flow over a two- 
dimensional body in an unbounded stream are: 

At infinity 


<J> = 0 

y 

and on the body 

V*< Y > - (2) 

It is now assumed that <t may be expanded in the form 

4> = x + 0^ + + . . . 

where 0 n+ q and its derivatives are small compared with 0 n , 0 n 

is of the order t n . Then, when the expansion for G> is substituted 
in the differential equation (1), the linear equations for 0-^ and 02 

are 

P 2 S^lxx + 01yy = 0 ( 3 ) 

P^02xx + 02yy = (l + Q )P^0 lyy0l y + 01xy01y] (*0 

O 

y + ± M^, p P 

where a = — P = 1 - Moo , Moo is the Mach number 1 at infinity, 

. 3 

and 7 is the ratio of specific heats at constant pressure and constant 
volume . 

The boundary conditions for 0q and 02 > to the order t and t 2 , 
respectively, are: 

On the body 

. , dY l 

01y(x,°) = t ^ ' (5) 

J o dY ? dY r 

02y(x,O) = t 2 + t 0 lx (x,O) + P^Y^^^O) (6) 

^lx “ ^ly = ^2x = ^2y = ® 


and at infinity 


( 7 ) 
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The well-known source-sink solution of the differential equation 
for valid for symmetric profiles and satisfying the boundary con- 

ditions given by equations (5) and ( 7 ), is (see reference 3, for example) 


0lU,y) 



logg^x - 4) 2 + 3 2 y 2 d£ 


( 8 ) 


The solution of the nonhomogeneous differential equation for 02, 
equation (4), can be expressed as the sum of the particular inte- 
gral ^(x^y) and a function <P2^ x ^y) satisfying the homogeneous 
equation. 


?2? 2xx + 


9 


2yy 


= o 


The particular integral ^ equation (4) (reference 2) is 


* 2 U,y) = M oo 2 0 ix 




(9) 


With the use of equations (8) and (9)> *2x and $ 2 y 8X6 easily shown 

to vanish at infinity and the boundary conditions for qp 2 , from equa- 
tions (6) and (7)> become: 


On the body 


9 2x (x,0) 


dF(x) 

dx 


where 


dF(x) 

dx 


and at infinity 


P 2t [^L0ixx( x > 0 ) + ^ 01x( x >°)] + t2 “ 
M tt 2 (l + |)0 lxy (x,O)0 1 (x,O) 


9 


2x 


- ”2y * 0 


It is to be understood that all differentiations axe made prior to 
setting y equal to zero. 
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The differential equation and boundary conditions for 92 , axe of 
the same form as for Therefore, the solution for 02 niay 

expressed as 


0 2 U,y) 


jit f di^ los e^ x ~ + p2y2 + V x ' y ) 


provided that 



|f( 6)| di converges and F(£) 


is sectionally con- 


tinuous. Since the source-sink solution may be obtained by using the 
concepts of Fourier integrals, it follows that the function F(£) must 
satisfy the conditions which are required for a function to be repre- 
sented in a Fourier integral. These conditions are sufficient to insure 
the existence of the solution although somewhat broader conditions 
for F(!) are known. These conditions seem to exclude the representa- 
tion of the second-order flow past blunt bodies by a continuous distri- 
bution of sources and sinks on the axis alone. 


Then the velocity potential to the order 


is 


®(x>y) 



dYi 

IT + 


dF(J) 

“dl 


l°ge\J(x - I ) 2 + d| 


+ 


^ g (x,y) + 0(t3j 


( 10 ) 


and the velocity components on the body to the same order are 
u = <D x (x,Y) 

= 1 + 


i P r 1 

. dYi 

df(l) 

d! 

J.x 

t IT + 

dl 

x - 1 


-TT + tYi0 lxy (x,O) + 




+ 0 (t 3 ) 


(ID 


v = «y(xyY) 

= t irf + ^ix (x >°l| +t2 ^r + 0(t3) (12) 

where P denotes the Cauchy principal value of the integral. 
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Once is determined for a given profile, the second-order 

velocity potential and the corresponding velocity components can be 
found from equations (10), (11), and (12). 

The third- and higher-order terms for the velocity potential or 
velocity components can be determined in the same manner as the second- 
order terms. The third-order term for the velocity potential is written 
as the sum of a function satisfying the homogeneous differential equa- 
tion and a particular integral (from reference 2). Then the boundary 
conditions for the homogeneous equation axe again in a form similar to 
those for 0i and the third-order solution is the sum of the particular 
integral and an expression for a source distribution with the source 
strength determined by the boundary conditions. 

It appears that this method, as developed herein, is not applicable 
for evaluating the flow over all bodies. For blunt bodies, the method 
must be refined to include other .solutions in conjunction with the 
source-sink solution in order to satisfy the boundary conditions * It 
might be noted that the troublesome terms in F(i) arise from the 
particular integral. Since the particular integral is multiplied by M 2 , 
the difficulty encountered is a Mach number effect. 

It does not seem possible to obtain a second-order solution for the 
flow past a sharp-nose body at a small angle of attack (to the second order 
in thickness and first order in angle of attack) by only a source-sink and 
vortex distribution on the axis. An integral equation must be solved to 
obtain the function in the lift problem corresponding to F( |) in the 
thickness problem. Therefore, it is difficult to make any general state- 
ments about the conditions this function must satisfy. However, the 
particular integral seems to introduce a singularity in the expression 
for the lift coefficient which cannot be canceled by the solution to the 
'homogeneous equation. Here again the method must be refined in order to 
obtain the second-order lift solutions. 

In general, the integral for 02 or its derivatives are more 
difficult to evaluate than those for 0q and in some cases it may be 
desirable to approximate the integral' for 02- When 01 is known, the 
integrand of the integral for 02 or its derivatives may be expanded 
in a Taylor series and the integration performed term t>y term. 

The Taylor expansion of the function ^ ) about the point x is 

d£ 

sli) .f - p( n+ i) ( , 

48 TSo n: 


( 13 ) 
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where 


- t2 sr + e 2 «i?W*,o) * u 2 t gi 0 lx (x,o) - 

M °° 2 ( 1 + |)^lxy( X ^ O )0i( X ^ 0 ) 


F (n) .(x) = F(x) 
dx 


Then, when - of equation (13) is substituted in equation (ll) and 

the resulting equation is integrated term by term, the u component of 
velocity becomes 


u = 1 + 0ixU,O) 


A 



-L dF(x) 
Pit dx 


logg 


1 - X 
1 + X 



(-1 < X < l) 


CALCULATION OF THE FLOW PAST THE KAPLAN BUMP 


As an example of the use of the second-order equations, the velocity 
at the surface of the Kaplan bump is evaluated. The equation of the 
bump is 

Y(x) = t(l - x 2 ) 3//2 - 3x 2 t 2 (l - x 2 ) 3/2 + 0(t3) 


Yq = (1 - x 2 ) 3 / 2 
Y 2 = -3x 2 (1 - X 2 ) 3 / 2 
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From equation (8), the expression for 0 -j y (x,O) is 

»1 


i (x0 ) = ipr^A = - 31 pr 

Hx ’ Pit J x - 5 Pit J _i 


s\H - r 


T 


d£ 


By the substitution Of I = cos 9, the expression for 0 lx (x,O) becomes 

**<*.<» -g(f ^ -.--r?*?,' 


p« 


X - cos 0, 


With sin 2 0 replaced by jj(l - cos 20) and with the use of the 
following integral 


cos n0 d0 -it sin na 


cos a - cos 0 sin a 


the expressions for 0 1y (x^O) and 0-^(x,O) are 

0 lx (x,°) = ^(| - x 2 ) 


0lU,O) = |(| x - x^J 

The u component of velocity (equation (ll)) is 


u = 4> x (x,Y) 


*’x' 

= 1 ' + 


p« ^ 


df(l) 


dS 


d| 


f + tY Axy (x ' 0) + 


(14) 


( -1 < x < 1) 


( -1 < x < 1) 


M 00 2 (l + |)|0 1 xx (x,O)0 1 (x,O) + ^(x^O)] U 0(t 3 ) 


( 15 ) 


dF(x) 

dx 


W>\ + t2 W- - 

M ro 2 (l + |)0 lxy (x,O)0 1 (x,O) 


with 
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d 2 Y 1 (x) 


where = t — — . The various terms arising in 

• dx^ 


dF(x) 


dx 


are 


6tx 


Y 101xx( x >°) = - p “ ' - p 


(1 - x 2 ) 3 / 2 = I^(x3 - x)\^-^ : 


dY- 

dx 


- 0ix(x,°) = -y(| - x 3 ) Vl - X 2 


a 

dx 


= (-6x + I5x 3 )\/i - 


,2 


d 2 Yi • 


dx 


- 0 1 (x,O) 


-3t^g x - 4x3 + 2x^j 


Jl - x 2 


y 


(16) 


The following integrals will be required: 

a 


P i t £ ZEm , P 




x - | 


0 


n cos 9 sin 2 0 d0 
x - cos 0 


m ”(-5 + x2 ) 


A 




■£ 


I dl 


= P 


= P 


' n cos 3 0 sin 2 0 d0 
, x - cos 9 


, 1 x 4 ’ 
rt ,_ 5 - y+ x 


cos 0 d0 


\[ r^~f(x - 1 ) Jo x • cos 0 


= -it 


'! *3 


l J dl 


= P 


,« 3 

rrtC'-J 


COS D Q d.0 


.1 - | 2 (x - |) J o x - cos 0 


■ -”(l + x2 ) 


■Ac: 


i 5 « 


„ r cos 0 d0 ( 3 x V 

= P / x = ■« ff + -s- + x 

Jq x - cos 0 18 2 


) 


rUT) 


ii - ? 2 (x - i) 

where repeated use has been made of equation (l4). 

From equation (15), with the use of equations (l6) and (17), wnH 


with 
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[^u,0j] 2 + ^(x,O)0 lxx (x,o) = pi(| - 6x 2 + 5X 11 ) 
the x-component of velocity on the body/ to the order t 2 , is 


u(x,Y) = 1 - | |(2 x 2 - l) + 2|_ 3X 1 * - 3x 2 + i + 5SX* 1 - § Ox 2 + / - 


I P 2 + M x*» - 2 x 2 + i 

8 00 '\2 2 4 , 

in agreement with reference 4 . 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., October 23, 1950 
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